THE INTERNAL DYNAMICS OF GLOBULAR CLUSTERS 



G. MEYLAN 

European Southern Observatory, Karl-Schwarzschild-Strasse 2, 
D-85748 Garching bei Miinchen, Germany 
E-mail: gmeylan@eso.org 

Galactic globular clusters are ancient building blocks of our Galaxy. They represent 
a very interesting family of stellar systems in which some fundamental dynamical 
processes have been taking place for more than 10 Gyr, but on time scales shorter 
than the age of the universe. In contrast with galaxies, these star clusters represent 
unique laboratories for learning about two-body relaxation, mass segregation from 
equipartition of energy, stellar collisions, stellar mergers, core collapse, and tidal 
disruption. This review briefly summarizes some of the tremendous developments 
that have taken place during the last two decades. It ends with some recent results 
on tidal tails around galactic globular clusters and on a very massive globular 
cluster in M31. 



1 Introduction 

There are about 150 globulars orbiting in the halo of our Galaxy. They look 
like huge swarms of stars, characterized by symmetry and apparent smooth- 
ness. Fig. 1 below displays an image of NGC 5139 = u Centauri, the brightest 
and most massive galactic globular cluster. This 40' by 40' image from the 
Digital Sky Survey does not reach, in spite of its rather large angular size, the 
outer parts of the cluster. With its tidal radius of about 40-50', the apparent 
diameter of uj Centauri on the plane of the sky is significantly larger than the 
apparent 30' diameter of the full moon. 

Globular clusters are old stellar systems, made of one single generation 
of stars. Although still somewhat uncertain, their individual ages range be- 
tween about 10 and 15 Gyr, with possible significant differences, up to a few 
gigayears, from one cluster to the other. Other properties of globular clusters 
exhibit significant variations: e.g., their integrated absolute magnitudes range 
from My 1 = -1.7 to -10.1 mag; their total masses from M to t = 10 3 to 5 x 
10 6 M Q ; their galactocentric distances from 2 to 120 kpc. 



2 A few dynamical time scales 

The dynamics of any stellar system may be characterized by the following 
three dynamical time scales: (i) the crossing time t cr , which is the time needed 
by a star to move across the system; (ii) the relaxation time t r i x , which is the 



meylanreview: submitted to World Scientific on February 1, 2008 1 



Figure 1. NGC 5139 = wCcntauri is the brightest and most massive galactic globular 
cluster. This image, from the Digital Sky Survey, has 40' by 40' in size. North is to the 
top, East to the left. 



time needed by the stellar encounters to redistribute energies, setting up a 
near-maxwellian velocity distribution; (hi) the evolution time t ev , which is 
the time during which energy-changing mechanisms operate, stars escape, 
while the size and profile of the system change. 

In the case of globular clusters, t cr ~ 10 6 yr, t r i x ~ 100 10 6 yr, and t ev <~ 
10 10 9 yr. It is worth mentioning that several (different and precise) defini- 
tions exist for the relaxation time. The most commonly used is the half-mass 
relaxation time t r h of Spitzcr (1987, Eq. 2-62), where the values for the mass- 
weighted mean square velocity of the stars and the mass density are those 
evaluated at the half-mass radius of the system (see Meylan & Heggie 1997 
for a review). It has been suggested that the combination of relaxation with 
the chaotic nature of stellar orbits in non-integrable potentials (e.g., most ax- 
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isymmetric potentials) causes a great enhancement in the rate of relaxation 
(Pfcnnigcr 1986, Kandrup & Willmcs 1994). Another suggestion which, if con- 
firmed, would revolutionise the theory of relaxation was made by Gurzadyan 
& Savvidy (1984, 1986) who proposed a much faster relaxation time scale 
than in standard theory, by a factor of order TV 2 / 3 . There is some support for 
this view on observational grounds (Vesperini 1992a, b). 

From size, luminosity, and mass points of view, globular clusters are brack- 
eted by open clusters on the lower side and dwarf elliptical galaxies on the 
upper side. Table 1 displays, for open clusters, globular clusters, and galax- 
ies, some interesting relations between the above three time scales. For open 
clusters, crossing time t cr and relaxation time t r i x are more or less equivalent, 
both being significantly smaller than the evolution time t ev . This means that 
most open clusters dissolve within a few gigayears. For galaxies, relaxation 
time t r i x and evolution time t ev are more or less equivalent, both being sig- 
nificantly larger than the crossing time t cr . This means that galaxies are not 
relaxed, i.e., not dynamically evolved. It is only for globular clusters that all 
three time scales are significantly different, implying plenty of time for a sig- 
nificant dynamical evolution in these stellar systems, although avoiding quick 
evaporation. 

Consequently, globular clusters represent an interesting class of dynamical 
stellar systems in which some dynamical processes take place on time scales 
shorter than their age, i.e., shorter than the Hubble time, providing us with 
unique dynamical laboratories for learning about two-body relaxation, mass 
segregation from equipartition of energy, stellar collisions, stellar mergers, and 
core collapse. All these dynamical phenomena are related to the internal dy- 
namical evolution only, and would also happen in isolated globular clusters. 
The external dynamical disturbances — tidal stripping by the galactic gravi- 
tational field — influence equally strongly the dynamical evolution of globular 
clusters. 
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3 Model building for globular clusters 

Already before the pioneering work of von Hocrncr (1960), who made the 
first iV-body calculations with N = 16, it was realized that computation of 
individual stellar motions could be replaced by statistical methods. Some 
parallels were drawn between a molecular gas and star clusters: the stars 
were considered as mass points representing the molecules in a collisionlcss 
gas. The analogy between a gas of molecules and a gas of stars is subject 
to criticisms, since the mean free path of a molecule is generally quite small 
compared with the size or scale height of the system, whereas the mean free 
path of a star is much larger then the diameter of the cluster; in addition 
molecules travel along straight lines, while stars move along orbits in the 
gravitational potential of all the other stars of the stellar system. Stellar 
collisions in clusters were studied by Jeans (1913), who remarked that they 
might be important in such stellar systems. The problem was then to seek 
the possible spherical distribution of such a gas in a steady state. 

3.1 Boltzmann's equation 

The commonest way of defining a model of a star cluster is in terms of its distri- 
bution function /(r, v, to), which is defined by the statement that fd 3 rd 3 vdm 
is the mean number of stars with positions in a small box d 3 r in space, ve- 
locities in a small box d 3 v and masses in an interval dm. In terms of this de- 
scription a fairly general equation for the dynamical evolution is Boltzmann's 
equation, 

|+v.V r /-V r$ .V v /=| ejic , (1) 

where $ is the smoothed gravitational potential per unit mass, and the right- 
hand side describes the effect of two-body encounters. The distribution / is 
a function of 7 variables if we take into account time. This is rather more 
than can be handle. But it is possible to reduce the complexity posed by 
Boltzmann's equation by taking moments. 

By taking moments of the Boltzmann's equation with respect to veloc- 
ities we obtain, for n = and 1, the Jeans equations which are expressions 
describing the rotation and the velocity dispersion: 

J Boltzmann ■ v™d 3 \r = JeansEqu. (2) 

By taking moments of the Jeans equations with respect to positions, we 
obtain the Tensor Virial equations which are expressions relating the global 
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kinematics to the morphology of the system, e.g., the ratio v /<r of ordered 
to random motions: 



In these ways, we obtain information about the general properties of so- 
lutions of Boltzmann's equation without recovering any solutions. 

3.2 Liouville's equation and Jean's theorem 

The general Boltzmann's equation can be greatly simplified in other ways. 
Because t cr is so short, after a few orbits the stars are mixed into a nearly 
stationary distribution, and so the term df/dt is practically equal to zero. 
In a similar way, because t r h is so long, the collision term (df/dt) enc can be 
ignored. What is left, i.e., 



is an equilibrium form of what is frequently called Liouville's equation, or the 
collisionless Boltzmann's equation, or the Vlasov equation. 

In simple cases, the general solution of Equ. 4 is given by Jeans' theorem, 
which states that / must be a function of the constants of the equations of 
motion of a star, e.g., of the stellar energy per unit mass e = v 2 /2 + $. 
Such quantities are also called integrals of the motion. If not all integrals 
of the motion are known, such functions are still solutions, though not the 
most general. For a self-consistent solution, the distribution function / must 
correspond to the density p required to provide the cluster potential <E> C , i.e.: 



Many different kinds of models may be constructed with this approach. In 
the first place there is considerable freedom of choice over which integrals of 
the motion to include. In the second place one is free to choose the func- 
tional dependence of these integrals, i.e., the analytic form of the distribution 
function (see, e.g., Binney 1982, and Binney & Tremaine 1987). 

King (1966) provided the first grid of models (with different concentra- 
tions c = log (r t /r c ) where r t and r c are the tidal and core radii, respectively) 
that incorporate the three most important elements governing globular cluster 
structure: dynamical equilibrium, two-body relaxation, and tidal truncation. 
These models depend on one integral of the motion only — the stellar energy 
per unit mass e — and the functional dependence is based on the lowered 
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(3) 



V.Vr/ - V r $.Vv/ = 



(4) 





maxwellian (see Equ. 6 below). Such models are spherical and their velocity 
dispersion tensor is everywhere isotropic. 

Models more complicated have been built since then. Da Costa & Free- 
man (1976) generalized the simple single-mass King models to produce more 
realistic multi-mass models with full equipartition of energy in the centre. 
Gunn & Griffin (1979) developed multi-mass models whose distribution func- 
tions depend on the stellar energy per unit mass e and the specific angular 
momentum I. Such models are spherical and have a radial anisotropic ve- 
locity dispersion (y% =/= v 2 . = v^). Called King-Michie models, they associate 
the lowered maxwellian of the King model with the anisotropy factor of the 
Eddington models: 

f(s, 1) <x (exp(-2j 2 £ ) - cxp(-2j 2 e t )) exp(-.7 2 Z 2 /r 2 ) (6) 

Lupton & Gunn (1987) developed multi-mass models whose distribution func- 
tions depend on a third integral of motion I 3 , in addition to the stellar energy 
per unit mass e and the component of angular momentum parallel to the 
rotation axis l z . Although no general analytical form for a third integral is 
available, the existence of an analytic third integral of motion I 3 in special 
cases has been known for decades, since the work by Jeans (1915). Because 
the rotation creates a non-spherical potential, I 3 = l 2 is in fact only an ap- 
proximate integral and Lupton & Gunn's distribution function does not obey 
the collisionless Boltzmann's equation for equilibrium (Eq. 4). 

These were notable landmarks in these developments, among many others. 
Table 2 hereafter, from Meylan & Hcggie (1997), list for the static models 
(King, King-Michie, 3-Integral) and for the evolutionary models (gas, Fokker- 
Planck, N-Body) the dynamical features and dynamical processes they take 
into account. Under the heading Dynamical Process, the second column in 
Table 2 states what kind of physical process it is that is named in the first 
column. 

4 Parametric and non-parametric approaches 

The method in the above section for analyzing globular cluster data is a model- 
building, or parametric, approach. One begins by postulating a functional 
form for the distribution function / and the gravitational potential $; often 
the two are linked via Poisson's equation, i.e. the stars described by / arc 
assumed to contain all of the mass that contributes to $. This / is then 
projected into observable space and its predictions compared with the data. 
If the discrepancies are significant, the model is rejected and another one is 
tried. If no combination of functions {/, <&} from the adopted family can be 
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Tabic 2. Dynamical models of globular star clusters 
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found that reproduces the data, one typically adds extra degrees of freedom 
until the fit is satisfactory. For instance, / may be allowed to depend on a 
larger number of integrals of the motion (Lupton & Gunn 1987) or the range 
of possible potentials may be increased by postulating additional populations 
of unseen stars (Da Costa & Freeman 1976). 
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This approach has enjoyed considerable popularity, in part because it is 
computationally straightforward but also because, as King (1981) has empha- 
sized, globular cluster data are generally well fitted by these standard models. 
But one never knows which of the assumptions underlying the models are 
adhered to by the real system and which are not. For instance, a deviation 
between the surface density profile of a globular cluster and the profile pre- 
dicted by an isotropic model is sometimes taken as evidence that the real 
cluster is anisotropic. But it is equally possible that the adopted form for 
/(e) is simply in error, since by adjusting the dependence of / on e one can 
reproduce any density profile without anisotropy. Even including the addi- 
tional constraint of a measured velocity dispersion profile does not greatly 
improve matters since it is always possible to trade off the mass distribution 
with the velocity anisotropy in such a way as to leave the observed dispersions 
unchanged (Dejonghe & Merritt 1992). Conclusions drawn from the model- 
building studies are hence very difficult to interpret; they are valid only to 
the extent that the assumed functional forms for / and $ are correct. 

These arguments suggest that it might be profitable to interpret kine- 
matical data from globular clusters in an entirely different manner, placing 
much stronger demands on the data and making fewer ad hoc assumptions 
about / and <f>. Ideally, the unknown functions should be generated non- 
parametrically from the data. Such an approach pioneered by Merritt (see, 
e.g., Merritt 1993a, b, 1996) has rarely been tried in the past because of the 
inherent instability of the deprojection process. We provide here after the 
results of two studies (parametric and non-parametric, respectively) of the 
globular cluster u> Centauri, both studies using exactly the same observational 
data (surface brightness profile and stellar radial velocities). 

4-1 Parametric approach applied to uj Centauri 

The mean radial velocities obtained with CORAVEL (Mayor et al. 1997) for 
469 individual stars located in the galactic globular cluster ui Centauri provide 
the velocity dispersion profile. It increases significantly from the outer parts 
inwards: the 16 outermost stars, located between 19.2' and 22.4' from the 
center, have a velocity dispersion a = 5.1 ± 1.6 kms -1 , while the 16 innermost 
stars, located within 1' from the center, have a velocity dispersion a = 21.9 ± 
3.9 kms -1 . This inner value of about er = 22 kins" 1 is the largest velocity 
dispersion value obtained in the core of any galactic globular cluster (Meylan 
etal. 1995). 

A simultaneous fit of these radial velocities and of the surface brightness 
profile to a multi-mass King-Michie dynamical model provides mean estimates 
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of the total mass equal to M to t = 5.1 10 6 Mq, with a corresponding mean 
mass-to-light ratio M/Ly =4.1. The present results emphasize the fact that 
u> Centauri is not only the brightest but also, by far, the most massive galactic 
globular cluster (Meylan etal. 1995). 

The fact that only models with strong anisotropy of the velocity dispersion 
(r a = 2-3 r c ) agree with the observations does not give a definitive proof 
of the presence of such anisotropy because of fundamental indetermination 
in the comparison between King-Michie models and observations. A strong 
anisotropy is nevertheless expected outside of the core of u> Centauri, given 
the large value of the half- mass relaxation time of about 26 < t r h< 46 10 9 yr) 
(Meylan etal. 1995). 

The reliability of the present application of King-Michie models might 
be questionable on a few fundamental points. In addition to the arbitrary 
choices of the two integrals of the motion and of the functional dependence 
of the distribution function on these two integrals, there is also the assump- 
tion of thermal equilibrium among the different mass classes in the central 
parts of the cluster. From a theoretical point of view, mass segregation has 
been one of the early important results to emanate from small N-body sim- 
ulations. Since then, large N-body simulations and models integrating the 
Fokker-Planck equation for many thousands of stars have fully confirmed the 
presence of equipartition. Thanks to the high angular resolution of the Hubble 
Space Telescope (HST) cameras (FOC and WFPC2), mass segregation has 
now been observed in the core of a few galactic globular clusters (see, e.g., 
Anderson 1997, 1999). In the case of 47Tucanae, the observed luminosity 
function by Anderson (1997, 1999) is in close agreement with equipartition- 
assuming King-Michie models and fails to fit the no-segregation models. This 
dichotomy is not as clear in the case of u> Centauri, probably because of its 
rather long central relaxation time. 

The problem about mass segregation does not concern its existence - 
it is happening — , but rather its quantitative evolution. Can there be an 
end to mass segregation, i.e., does the system ever reach a stable thermal 
equilibrium ? Underlying is the problem of core collapse (see Spitzer 1969, 
Chernoff & Weinberg 1990), which is briefly described in § 7 below. 

4-2 Non-parametric approach applied to ui Centauri 

The stellar dynamics of to Centauri is inferred from the same radial velocities 
of 469 stars used in § 4.1 (Mayor etal. 1997). By assuming that the resid- 
ual velocities are isotropic in the meridional plane, <r ro = a z = a, Merritt 
etal. (1997) derived the dependence of the two independent velocity disper- 
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Figure 2. Surface brightness (a) and space density (b) profiles for u> Centauri from Merritt 
ctal. (1997). 



sions a and on various positions in the meridional plane. The central 
velocity dispersion parallel to the meridional plane is cr — 171^6 kms -1 . 
With this approach, there is no evidence for significant anisotropy anywhere 
in uo Centauri. Thus, this cluster can reasonably be described as an isotropic 
oblate rotator (Merritt etal. 1997). 

The binned surface brightness measurements from Meylan (1986) are plot- 
ted in Fig. 2a, where the solid line in is an estimate of the surface brightness 
profile as the solution to the optimization problem. The estimate of 

space density profile v{r) may be defined as the Abel inversion of the estimate 



1 f c 
u{r)=— \ 

7T Jr 



dY, 



dR 



dR Vi? 2 - r 2 



(7) 



The dashed lines in Fig. 2b are 95% confidence bands on the estimate of v{r). 
Here r is an azimuthally-averaged mean radius. Both profiles are normalized 
to unit total number. This profile actually has a power-law cusp, v ~ r _1 , 
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inside of 0.5'; however the confidence bands are consistent with a wide range of 
slopes in this region, including even a profile that declines toward the center. 

The gravitational potential and mass distribution in lo Centauri are con- 
sistent with the predictions of a model in which the mass is distributed in the 
same way as the bright stars, although the cluster is assumed to be oblate 
and edge-on but mass is not assumed to follow light. The central mass den- 
sity is 2110j^gM©pc~ 3 . However this result may be strongly dependent on 
the assumption that the velocity ellipsoid is isotropic in the meridional plane. 
This central mass density determination is in full agreement with the values 
deduced from King-Michie models by Meylan etal. (1995). 

There is no significant evidence for a difference between the velocity dis- 
persions parallel and perpendicular to the meridional plane. The mass distri- 
bution inferred from the kinematics is slightly more extended than, though 
not strongly inconsistent with, the luminosity distribution. The derived two- 
integral distribution function f(e, l z ) for the stars in ui Centauri is fully con- 
sistent with the available data. 

Large amount of kinematical data (radial velocities and proper motions for 
a few thousand stars) will soon allow the efficient use of the non-parametric ap- 
proach in the case of the largest two galactic globular clusters, viz. lu Centauri 
and 47Tucanae (Freeman etal., Meylan etal., both in preparation). 

5 Systemic rotation of lo Centauri 

Systemic rotation in globular clusters has been expected for a long time, 
especially in lo Centauri, because of its significant flattening. The first clear 
evidence of such rotation was observed, in this cluster and in 47Tucanae, by 
Meylan & Mayor (1986). More recently, rather than fitting the data to a 
family of models, estimate of the rotation was obtain non-parametrically, by 
direct operation on the data by Merritt etal. (f997). 

Fig. 3 displays the contours of constant v<f, which are remarkably similar 
in shape to those of the parametric model postulated by Meylan & Mayor 
(1986), at least in the region near the center where the solution is strongly 
constrained by the data. The rotational velocity field is clearly not cylindrical; 
instead, has a peak value of 8 kms" 1 at about 7' from the center in the 
equatorial plane, and falls off both with increasing w and z. In the region 
inside the peak, the rotation is approximately solid-body; at large radii, the 
available data do not strongly constrain the form of the rotational velocity 
field. The mean motions are consistent with axisymmetry , once a correction 
is made for perspective rotation resulting from the cluster proper motion. 

The above inferred rotational velocity field in lo Centauri agree remark- 
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ably with the predictions of the theoretical model by Einsel & Spurzem (1999), 
who have investigated the influence of rotation on the dynamical evolution of 
collisional stellar systems by solving the orbit-averaged Fokker-Planck equa- 
tion in (e, / z )-space. However it is not clear that any relevant comparison 
exists: because of the long relaxation time in u> Centauri, the rotation prob- 
ably still reflects to a large extent the state of the cluster shortly after its 
formation. The observed estimate of v<f,(zu, z) might therefore be most useful 
as a constraint on cluster formation models. 

But things may be even more complicated ! Using their calcium abun- 
dances for about 400 stars with radial velocities by Mayor et al. (1997), Norris 
etal. (1997) found that the 20% metal-rich tail of the [Ca/H] distribution is 
not only more centrally concentrated, but is also kinematically cooler than 
the 80% metal-poor component. While the metal-poorer component exhibits 
well-defined systemic rotation, the metal-richer one shows no evidence of it, 
in contradistinction to the simple dissipative enrichment scenario of cluster 
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formation. 



6 Rotation vs. velocity dispersion 

All results about rotation depend on the value of the angle i between the plane 
of the sky and the axis of symmetry of the cluster. This angle remains un- 
known. Since the two best studied clusters, viz. u> Ccntauri and 47Tucanae, 
belong to the small group of clusters which, among the 150 galactic globular 
clusters, are the flattest ones, we can expect, from a statistical point of view, 
that their angles i should not be very different from 0° < i < 30°, the clus- 
ters being seen nearly edge-on. The importance of rotation (namely, of its 
projection along the line of sight) increases as i gets closer to 0°. 

The relative importance of rotational to random motions is given by the 
ratio v /(j 0l where vl is the mass- weighted mean square rotation velocity and 
(Tq is the mass- weighted mean square random velocity. For i = 90° and 60° , in 
lu Centauri the ratio v /a = 0.35 and 0.39 and in 47Tucanae the ratio v /a 
= 0.40 and 0.46, respectively (Meylan & Mayor 1986). Even with i = 45°, the 
dynamical importance of rotation remains weak compared to random motions. 
The ratio of rotational to random kinetic energies is ~ 0.1, confirming the fact 
that globular clusters are, above all, hot stellar systems. 

Rotation has been directly observed and measured in twelve globular clus- 
ters (see Table 7.2 in Meylan & Heggie 1997). The diagram (v /a vs. (e)), of 
the ratio of ordered v to random a Q motions as a function of the ellipticity (e) , 
has been frequently used for elliptical galaxies and its meaning is extensively 
discussed in Binney & Tremaine (1987 Chapter 4.3). The low luminosity (L 
^ 2.5 10 10 Lq) elliptical galaxies and spheroids have (v /<j 0l {e)) values which 
are scattered along the relation for oblate systems with isotropic velocity- 
dispersion tensors, while the high luminosity [L ^> 2.5 10 10 Lq) elliptical 
galaxies have (v /a ,(e)) values which are scattered below the above relation, 
indicating the presence of anisotropic velocity-dispersion tensors. Given their 
small mean ellipticities (0.00 < (e) < 0.12), globular clusters are located in 
the lower-left corner of the (v /a Q vs. (e)) diagram, an area characterized by 
isotropy or mild anisotropy of the velocity-dispersion tensor. 

7 Overwhole dynamical evolution towards core collapse 

Till the late ninety seventies, globular clusters were thought to be relatively 
static stellar systems since most surface-brightness profiles of globular clusters 
were successfully fitted by equilibrium models. Nevertheless, it had been 
already known, since the early sixties, that globular clusters had to evolve 
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dynamically, even when considering only relaxation, which causes stars to 
escape, consequently cluster cores to contract and envelopes to expand. But 
dynamical evolution of globular clusters was not yet a field of research by 
itself, since the very few theoretical investigations had led to a most puzzling 
paradox: core collapse (Henon 1961, Lynden-Bell & Wood 1968). 

It was only in the early eighties that the field grew dramatically. On 
the theoretical side, the development of high-speed computers allowed nu- 
merical simulations of dynamical evolution. Nowadays, Fokker-Planck and 
conducting-gas-sphere evolutionary models have been computed well into core 
collapse and beyond, leading to the discovery of possible post-collapse oscil- 
lations. In a similar way, hardware and software improvements of N-body 
codes provide very interesting first results for 10 4 -body simulations (Makino 
1996a, b, Spurzem & Aarseth 1996, Portegies Zwart etal. 1999), and give the 
first genuine hope, in a few years, for 10 5 -body simulations. On the obser- 
vational side, the manufacture of low-readout-noise Charge Coupled Devices 
(CCDs), combined since 1990 with the high spatial resolution of the Hubble 
Space Telescope (HST) , allow long integrations on faint astronomical targets 
in crowded fields, and provide improved data analyzed with sophisticated 
software packages. 

8 Gravothermal instability, gravothermal oscillations 

For many years (between about 1940 and 1960) secular evolution of globular 
cluster was understood in terms of the evaporative model of Ambartsum- 
ian (1938) and Spitzer (1940). In this model it is assumed that two-body 
relaxation attempts to set up a maxwellian distribution of velocities on the 
time scale of a relaxation time, but that stars with velocities above the es- 
cape velocity promptly escape. The next major step in understanding came 
when it was discovered that evolution arises also when stars escape from the 
inner parts of the cluster to larger radii, without necessarily escaping alto- 
gether. Antonov (1962) realised that these internal readjustments need not 
lead to a structure in thermal equilibrium, because thermal equilibrium may 
be unstable in self-gravitating systems (see Lynden-Bell & Wood 1968). The 
well known process of core collapse is interpreted as a manifestation of the 
gravothermal instability. 

Core collapse has been first observed and studied in simulations using 
gas and Fokker-Planck models. For an isolated cluster (without a tidal field) 
the time scale for the entire evolution of the core (when the density has for- 
mally become infinite) is about 15.7 t r h(0), when expressed in terms of the 
initial half-mass relaxation time (Cohn 1980). This result is for an isotropic 
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code starting from a Plummcr model with stars of equal mass, while for an 
anisotropic code the time extends to 17.6 t r h(0) (Takahashi 1995). 

The collapse time is generally shorter in the presence of unequal masses 
(Inagaki & Wiyanto 1984, Chernoff & Weinberg 1990). Murphy & Cohn 
(1988) give surface brightness and velocity dispersion profiles at various times 
during collapse, for a system with a reasonably realistic present-day mass 
spectrum. Addition of effects of stellar evolution, modeled as instantaneous 
mass loss at the end of main sequence evolution, delays the onset of core 
collapse (Angeletti & Giannone 1980, Applegate 1986, Chernoff & Weinberg 
1990, Kim etal. 1992). The effect of a galactic time-dependent tidal field can 
be to accelerate core collapse (Spitzer & Chevalier 1973). 

Examples of A-body models which illustrate various aspects of core col- 
lapse include Aarseth (1988), where TV = 1,000, Gicrsz & Heggie (1993) 
(N < 2, 000), Spurzem & Aarseth (1996) (N = 10, 000), and Makino (1996a,b; 
see Fig. 4 hereafter) (N < 32,000). 

At one time, it was not at all certain that a cluster could survive beyond 
the end of core collapse, with as a singularity characterized by infinite central 
density. Thus, many experts doubted whether the study of post-collapse 
clusters had any relevance to the interpretation of observations. Henon (1961, 
1965) showed that a cluster without such a singularity would evolve into one 
that did, and he realised that, in a real system, a flux of energy might well be 
supplied by the formation and evolution of binary stars, governing a series of 
collapses and expansions of the cluster core. 

Numerical simulations using gas and Fokker-Planck models show that sys- 
tems with at least a few thousand stars (Goodman 1987, Heggie & Ramamani 
1989, Breeden et al. 1994) follow a complicated succession of collapses and ex- 
pansions, called gravothermal oscillations by their discoverers (Sugimoto & 
Bettwieser 1983, Bettwieser & Sugimoto 1984). Quite apart from their rel- 
evance in nature, these oscillations are interesting in their own right, as an 
example of chaotic dynamics. From this point of view they have been studied 
by Allen & Heggie (1992), Breeden & Packard (1994), and Breeden & Cohn 
(1995). 

In 1995 the genuine occurrence of gravothermal oscillations in A-body 
systems was spectacularly demonstrated by Makino (1996a,b). These results 
confirm that the nature of post-collapse evolution in A-body systems is far 
more stochastic than in the simplified continuum models on which so much 
of our understanding rests at present. 
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Figure 4. Core collapse in systems with equal masses, from N-body integrations by Makino 
(1996b). The logarithm of the central density is plotted against time, scaled in proportion 
to the initial half-mass relaxation time. The successive curves, which correspond to different 
values of N, have been displaced vertically for clarity. For N < 32, 000 the first core collapse 
is clearly followed by gravothermal oscillations. 



9 Observational evidence of core collapse 

In the eighties, CCD observations allowed a systematic investigation of the 
inner surface brightness profiles (within ~ 3') of 127 galactic globular clus- 
ters (Djorgovski & King 1986, Chernoff & Djorgovski 1989, Trager etal. 
1995). These authors sorted the globular clusters into two different classes: (i) 
the King model clusters, whose surface brightness profiles resemble a single- 
component King model with a flat isothermal core and a steep envelope, and 
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(ii) the collapsed-core clusters, whose surface brightness profiles follow an al- 
most pure power law with an exponent of about -1. In the Galaxy, about 
20% of the globular clusters belong to the second type, exhibiting in their 
inner regions apparent departures from King-model profiles. Consequently, 
they are considered to have collapsed cores. 

The globular cluster M15 has long been considered as a prototype of the 
collapsed-core star clusters. High-resolution imaging of the centre of M15 
has resolved the luminosity cusp into essentially three bright stars. Post- 
refurbishment HST star-count data confirm that the 2.2" core radius observed 
by Lauer etal. (1991), and questioned by Yanny etal. (1994), is observed 
neither by Guhathakurta etal. (1996) with WFPC2 data nor by Sosin & 
King (1996) with FOC data. This surface-density profile clearly continues 
to climb steadily within 2". A maximum- likelihood method rules out a 2" 
core at the 95% confidence level. It is not possible to distinguish at present 
between a pure power-law profile and a very small core (Sosin & King 1996). 
Consequently, among the galactic globular clusters, M15 displays one of the 
best cases of clusters caught in a state of deep core collapse. 

10 Tidal tails from wide- field imaging 

10.1 Tidal truncation 

In addition to the effects of their internal dynamical evolution, globular clus- 
ters suffer strong dynamical evolution from the potential well of their host 
galaxy (Gncdin & Ostriker 1997, Murali & Weinberg 1997). These external 
forces speed up the internal dynamical evolution of these stellar systems, ac- 
celerating their destruction. Shocks are caused by the tidal field of the galaxy: 
interactions with the disk, the bulge and, somehow, with the giant molecular 
clouds, heat up the outer regions of each star cluster. The stars in the halo are 
stripped by the tidal field. All globular clusters are expected to have already 
lost an important fraction of their mass, deposited in the form of individual 
stars in the halo of the Galaxy (see Meylan & Heggie 1997 for a review). 

Recent N-body simulations of globular clusters embedded in a realistic 
galactic potential (Oh & Lin 1992; Johnston etal. 1999) were performed in 
order to study the amount of mass loss for different kinds of orbits and different 
kinds of clusters, along with the dynamics and the mass segregation in tidal 
tails. Grillmair etal. (1995) in an observational analysis of star counts in the 
outer parts of a few galactic globular clusters found extra-cluster overdensities 
that they associated partly with stars stripped into the Galaxy field. 
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10.2 Tidal tails from wide-field observations 

Leon, Meylan & Combes (2000) studied the 2-D structures of the tidal tails 
associated with 20 galactic globular clusters, obtained by using the wavelet 
transform to detect weak structures at large scale and filter the strong back- 
ground noise for the low galactic latitude clusters. They also present N-body 
simulations of globular clusters in orbits around the Galaxy, in order to study 
quantitatively and geometrically the tidal effects they encounter (Combes, 
Leon & Meylan 2000). 

Their sample clusters share different properties or locations in the Galaxy, 
with various masses and structural parameters, ft is of course necessary to 
have very wide field imaging observations. Consequently, they obtained, dur- 
ing the years 1996 and 1997, photographic films with the ESO Schmidt tele- 
scope. The field of view is 5.5° x 5.5° with a scale of 67.5"/mm. The filters 
used, viz. BG12 and RG630, correspond to B and R, respectively. All these 
photographic films were digitalized using the MAMA scanning machine of the 
Observatoire de Paris, which provides a pixel size of 10 fim. The astrometric 
performances of the machine are described in Bcrgcr ct al. (1991). 

The next step — identification of all point sources in these frames — was 
performed using SExtractor (Bertin & Arnouts 1996), a software dedicated to 
the automatic analysis of astronomical images using a multi-threshold algo- 
rithm allowing good object deblcnding. The detection of the stars was done 
at a 3-cr level above the background. This software, which can deal with huge 
amounts of data (up to 60,000 x 60,000 pixels) is not suited for very crowded 
fields like the centers of the globular clusters, which were simply ignored. A 
star/galaxy separation was performed by using the method of star/galaxy 
magnitude vs. log(star/galaxy area). 

For each field, a (B vs. B — V) color-magnitude diagram was constructed, 
on which a field/cluster star selection was performed, following the method 
of Grillmair etal. (1995), since cluster stars and field stars exhibit different 
colors. In this way present and past cluster members could be distinguished 
from the fore- and background field stars by identifying in the CMD the area 
occupied primarily by cluster stars. The envelope of this area is empirically 
chosen so as to optimize the ratio of cluster stars to field stars in the relatively 
sparsely populated outer regions of each cluster. 

10.3 Wavelet Analysis 

With the assumption that the data can be viewed as a sum of details with 
different typical scale lengths, the next step consists of disentangling these 
details using the space-scale analysis provided by the Wavelet Transform (WT, 
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cf. Slczak etal. 1994; Resnikoff & Wells 1998). Any observational signal 
includes also some noise, which has a short scale length. Consequently the 
noise is higher for the small scale wavelet coefficients. Monte-Carlo simulations 
were performed to estimate the noise at each scale and apply a 3-cr threshold 
on the wavelet coefficients to keep only the reliable structures. In this way it is 
possible to subtract the short- wavelength noise without removing details from 
the signal which has longer wavelengths. The remaining ovcrdensities of the 
cluster-like stars, remaining after the application of the wavelength transform 
analysis to the star counts, are associated with the stars evaporated from the 
clusters because of dynamical relaxation and/or tidal stripping by the galactic 
gravitational field. 

It is worth emphasizing that in this study, the following strong observa- 
tional biases were taken into account: (i) bias due to the clustering of galactic 
field stars; (ii) bias due to the clustering of background galaxies; (iii) bias due 
to the fluctuations of the dust extinction, as observed in the IRAS 100-/xm 
map. 

10. 4 Observational Results 

The most massive galactic globular cluster, wCentauri (Meylan etal. 1995), 
currently crossing the disk plane, is a nearby globular cluster located at a 
distance of 5.0 kpc from the sun. Its relative proximity allows, for the star 
count selection, to reach the main sequence significantly below the turn-off. 
Estimates, taking into account the possible presence of mass segregation in 
its outer parts, show that about 0.6 to 1 % of its mass has been lost during 
the current disk shocking event. Although this cluster has, in this study, one 
of the best tail/background S/N ratios, it is by far not the only one exhibiting 
tidal tails. 

Considering all 20 clusters of the sample, the following conclusions are 
reached (see Leon, Meylan & Combes 2000 for a complete description of this 
work): 

• All the clusters observed, which do not suffer from strong observational 
biases, present tidal tails, tracing their dynamical evolution in the Galaxy 
(evaporation, tidal shocking, tidal torquing, and bulge shocking). 

• The clusters in the following sub-sample (viz. NGC 104, NGC 288, 
NGC 2298, NGC 5139, NGC 5904, NGC 6535, and NGC 6809) exhibit 
tidal extensions resulting from a recent shock, i.e. tails aligned with the 
tidal field gradient. 

• The clusters in another sub-sample (viz. NGC 1261, NGC 1851, 
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Figure 5. NGC 5139 = ui Ccntauri. In the upper panel, filtered image of color-selected 
star-count overdensities using the Wavelet Transform to be compared with the raw star 
counts in the lower panel. The upper panel displays the full resolution using the whole set 
of wavelet planes. From Leon etal. (2000). 
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NGC 1904, NGC 5694, NGC 5824, NGC 6205, NGC 7492, Pal 5, and 
Pal 12) present extensions which are only tracing the orbital path of the 
cluster with various degrees of mass loss. 

• NGC 7492 is a striking case because of its very small extension and 
its high destruction rate driven by the galaxy as computed by Gnedin & 
Ostriker (1997). Its dynamical twin for such an evolution, namely Pal 12, 
exhibits, on the contrary, a large extension tracing its orbital path, with 
a possible shock which happened more than 350 Myr ago. 

• The presence of a break in the outer surface density profile is a reliable 
indicator of some recent gravitational shocks. 

Recent CCD observations with the Wide Field Imager at the ESO/MPI 
2.2-m telescope and with the CFH12K camera at the Canada-France-Hawaii 
3.6-m telescope will soon provide improved results, because of the more ac- 
curate CCD photometry. These observations will allow more precise observa- 
tional estimates of the mass loss rates for different regimes of galaxy-driven 
cluster evolution. 

10.5 Numerical Simulations 

Extensive numerical N-body simulations of globular clusters in orbit around 
the Galaxy were performed in order to study quantitatively and geometrically 
the tidal effects they encounter and to try to reproduce the above observations. 
The N-body code used is an FFT algorithm, using the method of James (1977) 
to avoid the periodic images. With N = 150,000 particles, it required 2.7 
seconds of CPU per time step on a Cray-C94. 

The globular clusters are represented by multi-mass King-Michie models, 
including mass segregation at initial conditions. The Galaxy is modeled as 
realistically as possible, with three components, bulge, disk and dark halo: 
the bulge is a spherical Plummer law, the disk is a Miyamoto-Nagai model, 
and the dark matter halo is added to obtain a flat galactic rotation curve. 

The main conclusions of these simulations can be summarized as follows 
(see Combes, Leon & Meylan 2000 for a complete description of this work): 

• All runs show that the clusters are always surrounded by tidal tails and 
debris. This is also true for those that suffered only a very slight mass loss. 
These unbound particles distribute in volumic density like a power-law as 
a function of radius, with a slope around -4. This slope is much steeper 
than in the observations where the background-foreground contamination 
dominates at very large scale. 
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Figure 6. Tidal tails mapped at different epochs with the wavelet algorithm applied to 
one of our simulations. The direction perpendicular to the galactic plane is indicated by 
the arrow. The time sequence starts with the lower-left panel and ends with the upper 
right one. The third panel exhibits tails which are quite reminiscent of what is observed in 
NGC 5139 = wCentauri (see Fig. 5 above). From Combes etal. (2000). 



• These tails arc preferentially composed of low mass stars, since they are 
coming from the external radii of the cluster; due to mass segregation 
built up by two-body relaxation, the external radii preferentially gather 
the low mass stars. 
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• For sufficiently high and rapid mass loss, the cluster takes a prolate shape, 
whose major axis precesses around the z-axis. 

• When the tidal tail is very long (high mass loss) it follows the cluster 
orbit: the observation of the tail geometry is thus a way to deduce cluster 
orbits. Stars are not distributed homogeneously through the tails, but 
form clumps, and the densest of them, located symmetrically in the tails, 
are the tracers of the strongest gravitational shocks. 

Finally, these N-body experiments help to understand the recent observations 
of extended tidal tails around globular clusters (Grillmair et al. 1995, Leon 
ct al. 2000) : the systematic observations of the geometry of these tails should 
provide much information on the orbit, dynamics, and mass loss history of 
the clusters, and on the galactic structure as well. 

11 Gl in M31: globular cluster or dwarf galaxy ? 

The globular cluster Mayall II = Gl, recently observed with the Hubble Space 
Telescope (HST) camera WFPC2 (Rich etal. 1996, Jablonka et al. 1999, 2000, 
Meylan etal. 2000), is a bright star cluster which belongs to our companion 
galaxy, Andromeda = M31. Its integrated visual magnitude V — 13.75 mag 
corresponds to an absolute visual magnitude My = -10.86 mag, with E(B— V) 
= 0.06 and a distance modulus (m — M)mzi = 24.43 mag, implying a total 
luminosity of about Ly ~ 2 x 10 6 L Q . 

The coordinates of Gl, viz. a G i(J2000.00) = 00° 32' 46.878" and 
6 G i(J2000.00) = +39° 34' 41.65", when compared to the coordinates of the 
center of M31, viz. a M 3i(J2000.00) = 00° 42' 44.541" and <5 M 3i (J2000.00) = 
+41° 16' 28.77", place it at a projected distance of about 3°, i.e. 39.5 kpc 
from the center of M31. In spite of this rather large projected distance, both 
color-magnitude diagrams and radial velocities of Gl and M31, viz. V r (Gl) 
= - 331 ± 24 kms" 1 while V r (M31) = - 300 ± 4 kms" 1 (21-cm HI line) and 
1^(M31) = - 295 ± 7 kins" 1 (optical lines), completely support the idea that 
this cluster belongs to the globular cluster system of M31. 

Our (V vs. V —I) color-magnitude diagram reaches stars with magnitudes 
fainter than V = 27 mag, with a well populated red horizontal branch at about 
V = 25.25 mag; we confirm the existence of a blueward extension of the red 
horizontal branch clump as already observed by Rich ct al. (1996). From 
model fitting, we determine a rather high mean metallicity of [Fe/H] = -0.95 
± 0.09, somewhat between the previous determinations of [Fe/H] = -0.7 (Rich 
et al. 1996) and [Fe/H] = -1.2 (Bonoli 1987; Brodie & Huchra 1990). 
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From artificial star experiments, in order to estimate our true measure- 
ment errors, we observe a clear spread in our photometry that we attribute 
to an intrinsic metallicity dispersion among the stars of Gl. Namely, adopt- 
ing E(V — I) = 0.10 implies a 1-cr [Fc/H] dispersion of ± 0.50 dex; adopting 
E(V — I) = 0.05 implies a l-er [Fc/H] dispersion of ± 0.39 dex. In all cases, 
the intrinsic metallicity dispersion is significant and may be the consequence 
of self enrichment during the early stellar/dynamical evolution phases of this 
cluster. 

We have at our disposal two essential observational constraints allowing 
the mass determination of Mayall II = Gl: 

(i) First, its surface brightness profile from HST/WFPC2 images, pro- 
viding essential structural parameters: the core radius r c = 0.14" = 0.52 pc, 
the half-mass radius = 3.7" = 14 pc, the tidal radius r< ~ 54" = 200 pc, 
implying a concentration c = log {r t /r c ) ~ 2.5 (Meylan etal. 2000). 

(ii) Second, its central velocity dispersion from KECK/HIRES spectra, 
providing an observed velocity dispersion a b s — 25.1 kms^ 1 , and an aperture- 
corrected central velocity dispersion a = 27.8 kms -1 . 

11.1 King model and Virial mass estimates 

We can first obtain simple mass estimates from King models and from the 
Virial (see, e.g., Illingworth 1976). The first estimate, King mass, is given by 
the simple equation: 

King mass = p c r c (i — 167 r c n,a1 (8) 

where the core radius r c = 0.52 pc, the dimensionless quantity fi = 220 for c 
= log (r t /r c ) = 2.5 (King 1966), and the central velocity dispersion a = 27.8 
kms -1 . These values determine a total mass for the cluster of M tot = 15 x 
10 6 M Q with the corresponding M/L v ~ 7.5. 

The second estimate, Virial mass, is given by the simple equation: 

Virial mass = 670 r h <r 2 (9) 

where the half-mass radius = 14 pc and central velocity dispersion <r = 
27.8 kms' 1 . These values determine a total mass for the cluster of M tot — 
7.3 x 10 6 M Q with the corresponding M/L v ~ 3.6. 

11.2 King-Michie model mass estimate 

The existing observational constraints allow the use of a multi-mass King- 
Michie model as defined by Equ. 6 above. See §4.1 above and Meylan etal. 
(1995) in the case of such a model applied to wCentauri. In the case of 
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Figure 7. Surface brightness profile of the globular cluster Mayall II = Gl, from 
HST/WFPC2 shallow and deep images in F555W ~ V filter; the continuous line repre- 
sents a King-Michic model (first model in Table 3) fitted to the observed profile (Meylan 
ctal. 2000). 
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Table 3. Multi-mass King-Michic models for Mayal II = Gl 
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Gl, such a model is simultaneously fitted to the surface brightness pro- 
file from HST/WFPC2 and to the central velocity dispersion value from 
KECK/HIRES. 

An extensive grid of about 150,000 models was computed in order to 
explore the parameter space defined by the Initial Mass Function (IMF) ex- 
ponent x, where x would equal 1.35 in the case of Salpeter (1955), the central 
gravitational potential W , and the anisotropy radius r a . The IMF exponent 
consists actually of three parameters, Xh r , describing the heavy remnants, re- 
sulting from the already evolved stars with initial masses in the range between 
0.85 and 100 M Q ; x^ s , describing the stars still on the Main Sequence, with 
initial masses in the range between 0.25 and 0.85 M Q ; and xf}g n describing 
the stars still on the Main Sequence, with initial masses in the range between 
0.10 and 0.25 M . 

Table 3 presents eleven of the 50 models with the lowest % 2 , illustrating 
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some of the input and output parameters. Good models are considered as such 
not only on the basis of the x 2 of the surface brightness fit (see Fig. 7) , but 
also from their predictions of the observed integrated luminosity of the cluster 
and of the input mass-to-light ratio of the model. The different columns in 
Table 3 give, for each model, its IMF exponents x^ s and xf}g n ; the fraction 
Mhr of its total mass in the form of heavy stellar remnants such as neutron 
stars and white dwarfs; its concentration c = log (r t /r c ); its total mass M to t 
of the cluster, in solar units; and its corresponding mass-to-light ratio M/Ly 
also in solar units. Since the velocity dispersion profile is reduced to one 
single value — the central velocity dispersion — the models are not strongly 
constrained, providing equally good fits to rather different sets of parameters. 

The IMF exponent Xh r , describing the amount of neutron stars, appears 
in all models to be very close to x = 1.35 (Salpeter 1955). Given the lack 
of constraint from the absence of any velocity dispersion profile, the most 
reliable results are related to the concentration and the total mass. With a 
concentration c = log (r t /r c ) somewhere between 2.45 and 2.65, Gl presents 
clearly and in all cases the characteristics of a collapsed cluster. This is com- 
pletely different from uj Centauri, the most massive but losse galactic globular 
cluster, which, with a concentration of about 1.3, has a very large core ra- 
dius of about 5 pc and is consequently very far from core collapse. With a 
total mass somewhere between 13 and 18 10 6 M Q , and with the corresponding 
mass-to-light ratio MjLy between 6 and 9, Gl is significantly more massive 
than u> Centauri, maybe by up to a factor of three. The King-Michie mass 
estimates are in full agreement with the King mass estimate, while the Virial 
mass estimate is smaller by about a factor of two. It is worth mentioning 
that such a mass difference is not typical of Gl: the same factor of about two 
is also observed between the King-Michie and Virial mass estimates of any 
cluster. See, e.g., Meylan & Mayor (1986) and Meylan et al. (1995) in the case 
of u> Centauri. 

11.3 Mayall 11= Gl is a genuine globular cluster 

From these three various mass determinations (King, Virial, King-Michie), we 
can reach the following conclusions about Mayall II = Gl: 

(i) All mass estimates give a total mass up to three times as large as the 
total mass of u> Centauri; 

(ii) With c = log (r t /r c ) ~ 2.5, Gl is more concentrated than 47Tucanae, 
which is a massive galactic globular cluster considered on the verge of collaps- 
ing; Gl has a surface brightness profile typical of a collapsed cluster; 

(iii) Gl is the heaviest of the weighted globular clusters. 
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Given these results we can wonder if, even more than wCentauri, Gl 
could be a kind of transition step between globular clusters and dwarf elliptical 
galaxies. There is a way of checking this hypothesis. Kormendy (1985) used 
the four following quantities — the central surface brightness ji , the central 
velocity dispersion a , the core radius r c , and the total absolute magnitude M 
— in order to define various planes from combinations of two of the above four 
quantities, e.g., (/i vs. log r c ). In all these planes, the various stellar systems 
plotted by Kormendy (1985) segregate into three well separated sequences: (i) 
ellipticals and bulges, (ii) dwarf ellipticals, and (iii) globular clusters. When 
plotted on any of these planes, Gl appears always on the sequence of globular 
clusters, and cannot be confused or assimilated with either ellipticals and 
bulges or dwarf ellipticals. The same is true for ui Centauri. 

Consequently, Mayall II = Gl can be considered a genuine bright and 
massive globular cluster. Actually, Gl may not be the only such massive 
globular cluster in M31. This galaxy, which has about twice as many globular 
clusters as our Galaxy, has at least three other clusters with central velocity 
dispersion larger than 20 kms -1 (Djorgovski etal. 1997). Unfortunately, so 
far, Gl is the only such cluster imaged with the high spatial resolution of the 
HST/WFPC2 camera, and consequently the only such massive cluster with 
known structural parameters. Gl and the other three bright M31 globular 
clusters represent probably the high-mass and high-luminosity tails of the 
otherwise very normal mass and luminosity distributions of the rich M31 
population of globular clusters. 

12 Conclusion 

This review summarizes only parts of the tremendous developments that have 
taken place during the last two decades. These recent developments are 
far from having exploited all the new capabilities offered by the impressive 
progress in computer simulations, made possible by more powerful single- 
purpose hardware and software (Hut & Makino 1999, Spurzem 1998). 

Observations too have still a lot of information to provide, which will 
require more elaborate modeling before full interpretation is reached. The 
mere observation of globular cluster stellar populations presents some puzzles 
which are far from being understood (Anderson 1997, 1999). 

The kinematical and dynamical understanding of globular clusters will 
need the exploitation of numerous radial velocities and proper motions of in- 
dividual stars. Only small quantities of radial velocities have been painfully 
accumulated over the last two decades, while the proper motions have so far 
simply been ignored. But there is an enormous amount of untapped infor- 
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mation locked in the radial velocities (for one third) and proper motions (for 
two thirds). Fortunately, a large amount of kinematical data (radial velocities 
and proper motions for a few thousand stars) will soon permit investigation 
of the 3-D space velocity distribution and rotation in the two largest galac- 
tic globular clusters, viz. u Ccntauri and 47Tucanae (Freeman etal., Meylan 
etal., both in preparation). 
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